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In the heavy-quark limit, the valence Fock-state components in the B mesons are de¬ 
scribed by a set of two light-cone wavefunctions. We show that these two wavefunctions 
obey simple coupled differential equations, which are based on the equations of motion 
in the Heavy Quark Effective Theory (HQET), and the analytic solutions for them are 
obtained. The results generalize the recently obtained longitudinal-momentum distribu¬ 
tion in the Wandzura-Wilczek approximation by including the transverse momenta. We 
find that the transverse momentum distribution depends on the longitudinal momentum 
of the constituents, and that the wavefunctions damp very slowly for large transverse 
separation between quark and antiquark. 

Keywords: B meson; light-cone wavefunctions; transverse momentum; heavy quark ef¬ 
fective theory. 

Along with the progress in both theory and experiment, B physics becomes one of 
the most active research areas in high energy physics. Many B meson exclusive de¬ 
cay processes turn out to be calculable systematically in the frameworks of newly de¬ 
veloped factorization formalisms, the so-called pQCD approach and QCD Factor¬ 
ization approach.In all of the calculations based on the factorization approaches, 
the light-cone distribution amplitudes of the participating mesons, which express the 
nonperturbative long-distance contributions in the factorized amplitudes,play an 
important role in making reliable predictions. The light-cone distribution ampli¬ 
tudes describe the probability amplitudes to find particular partons with definite 
light-cone momentum fraction in a meson, and thus are process-independent quan- 
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tity. It is well-known that, for the light mesons, the model-independent framework 
to construct the light-cone distribution amplitudes is established for leading and 
higher twists as well.®’® However, unfortunately, the distribution amplitudes for 
the B mesons are not well-known at present and they provide a major source of 
uncertainty in the calculations of the decay rates. 

Recently,^® we have presented the first systematic study for the B meson light- 
cone distribution amplitudes, and derived explicit forms for the quark-antiquark 
distribution amplitudes, which exactly satisfy the constraints coming from the equa¬ 
tions of motion and heavy-quark symmetry. We have found that the “Wandzura- 
Wilczek-type” contributions which correspond to the valence quark distributions, 
are determined uniquely in analytic form in terms of A, a fundamental mass param¬ 
eter of Heavy Quark Effective Theory (HQET).^^’^® We have also shown that both 
leading- and higher-twist distributions receive the contributions from the multi¬ 
particle states with additional dynamical gluons, and derived the exact integral 
representations for these contributions. 

By definition, the light-cone distribution amplitudes are given by the light-cone 
(Bethe-Salpeter) wavefunctions at zero transverse separation of the constituents. 
Thus the previous results of Ref.^® have been obtained for the configuration in 
which the quark and antiquark are separated by exactly light-like distance. The in¬ 
formation on transverse momentum distribution has been integrated out. However, 
the light-cone wavefunctions with transverse momentum dependence are necessary 
for computing the power corrections to the exclusive amplitudes, and also for esti¬ 
mating the transition form factors for B ^ D, B t:, etc, which constitute another 
type of long-distance contributions appearing in the factorization approaches for the 
exclusive B meson decays. 

In this Letter, we extend the analysis of Ref.^® to include the transverse mo¬ 
mentum effects. In particular, we derive explicit analytic formulae for a com¬ 
plete set of the B meson light-cone wavefunctions within the Wandzura-Wilczek 
approximation.^® It should be noted that the Wandzura-Wilczek approximation is 
not equivalent to the free field approximation. The leading Eock-states, which cor¬ 
respond to the twist-2 contribution in the case of light meson wavefunctions, carries 
the effect from the QCD interaction.®’® We also estimate the effects neglected in 
this “valence” approximation; we derive the exact result for the first moment of 
the transverse momentum squared in terms of the full light-cone wavefunctions, 
which include the higher Eock-states with additional dynamical, nonperturbative 
gluons.®' 

The light-cone wavefunctions are related to the usual Bethe-Salpeter wavefunc¬ 
tions at equal light-cone time z~^ = ( 2 ® -I- z®)/-\/2.^ The quark-antiquark light-cone 
wavefunctions '4’±{t,z^) of the B mesons in the heavy-quark limit are defined by 
the vacuum-to-meson matrix element of nonlocal operators, following Refs.^®’^^: 

{0\q{z)Thy{0)\B{p)) 

''^Radiative corrections due to the virtual gluons and/or quark-antiquark pairs can be included as 
the renormalization scale-depenence of the wavefunctions, which is governed by the renormalization 
group equations for the nonlocal operators in Eq. (1) below. The discussion of this point is beyond 
the scope of this work. 
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Here = 0, v'^ = 1, t = v ■ z, and = Mv^ is the 4-momentum of the B meson 
with mass M. hy(x) denotes the effective 6-quark field, h{x) « ex.p{—imtv-x)h^{x), 
and is subject to the on-shell constraint, fihy = F is a generic Dirac matrix 

and, here and in the following, the path-ordered gauge factors are implied in between 
the constituent fields, fs is the decay constant defined as usually as 


(0|(7(0)7^756„(0)|H(p)) =ifBMv^ , 


( 2 ) 


so that V'±(t = 0, = 0) = 1. Eq. (1) is the most general parameterization 

compatible with Lorentz invariance and the heavy-quark limit. 

Note that in the definition (1) the separation z^ = (0, z~, zt) between quark and 


antiquark is not restricted on the light-cone, z^ 




unlike the corresponding 


definition of the distribution amplitudes.Thus the distribution amplitudes 
4>±{t) in the notation of Refs.^°’^®’^® are given by the light-cone limit of the above 
wavefunctions as (j)±{t) = tl>±{t,z^) 

2^ —^0 

We also introduce the Fourier transforms with respect to the longitudinal sepa¬ 
ration t by 


'4’±{tjz'^)= j doj e z^) 


( 3 ) 


Here uiv'^ has the meaning of the light-cone projection fc'*" of the light-antiquark 
momentum in the B meson. 

We now exploit the constrains from the equations of motion. The procedure 
here is completely in parallel with that of our previous work, so we refer the readers 
to Ref.^° for the detail. The matrix elements of the exact operator identities (Eqs. 
(3), (4) of Ref.^°) from the equations of motion yield a system of four differential 
equations: 


d^p- 

/ dil!+ dip- 


(w-A) 


(w-A) 


d^P+ 

du 


dtp- 

92 dip+ 


,2 


+ — 0 ) 


+ 2 


9^2 9z2 
92 dlP+ 


+ 2{'ip+- = 0 , 


1 


90.2 5,2 +2(3V'+-V'-)=0, 


fdi^+ 

dp^-\ 

-1- 9 

(di^+ 

d-ip- \ 

1 doj 

doj J 


1 9z2 

9z2 ) 


( 4 ) 

( 5 ) 

( 6 ) 

0. (7) 


Here, 


A = M — mb 


iv ■ d{0\qThy\B{p)) 
{0\qrhy\Bip)) 


( 8 ) 


is the usual “effective mass” of meson states in the HQET,^2.i7 shorthand 

notation tp± = ■0±(a., z2) is adopted. These Eqs. (4), (5), (6), and (7) correspond 
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to Eqs. (7), (8), (10), and (11) in Ref.^°, respectively, but are given in the “w- 
representation” instead of the “<-representation” via Eq. (3). Further differences 
compared with Ref.^° are on two aspects. Firstly, the light-cone limit is not taken 
in order to explore the tranverse momentum distribution, so that the terms propor¬ 
tional to appear in the above Eq. (4). Secondly, we restrict our interest within 
the two-particle Fock states neglecting the contribution from the quark-antiquark- 
gluon three-particle operators. 

From Eqs. (5) and (6), we obtain 

^ -b V'+ + V-- = 0 . (9) 

This can be integrated with boundary conditions 'ip±(to, z^) = 0 for w < 0 or w ^ oo 
as 

(w — 2A) = 0 . (10) 

A system of Eqs. (4) and (10) has been solved for ^ 0 and completely determined 
the wavefunctions in the light-cone limit.The results are 

V’+(u;,z^ = 0) = (l>+{uj) = ^9{uj)9{2A - uj) , (11) 

V'-(w,z^ = 0) = = ^^^^9{uj)9{2A-uj) . (12) 

These solutions for z^ = 0 serve as “boundary conditions” to solve Eqs. (4)-(7) for 
z^ yb 0. From Eqs. (4) and (10) for z^ ^ 0, we find 

■0±(w,z^) = (()±(w)^(z^u;( 2A-w)) , (13) 

where ^{x) is some function of a single variable x, and satisfies ^(0) = 1 due to Eqs. 
(11), (12). The functional form of ^(x) can be easily determined from a remaining 
differential equation, e.g., (6), which was useless in the light-cone limit^° (Eq. (7) 
gives the result identical to Eq. (6)). We obtain ^(x) = Jq [y/—x), where Jq is a 
(regular) Bessel function, so that the analytic solution for the coupled differential 
equations (4)-(7) is given by 

%P±{uj,-Zt) = 4>±{^)Jo [\zt\\Iuj{ 2A - Uj'^ . (14) 

These are the light-cone wavefunctions for the transverse separation between 
quark and antiquark. 

For the momentum-space wavefunctions kx) defined by 

^±{t,-z^x) = J dwd^kr A;t) , (15) 

our solution gives 

^+{uj,kT) = ■^^^9{uj)9{2A-uj)6(k‘^-uj{2A-uj)) , (16) 

'il)-{oj,kT) = - j^9{uj) 9{2A — uj)5 (kx — u){ 2A — u))) . (17) 

27rA^ 
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The results (16) and (17) give exact description of the valence Fock components of 
the B meson wavefunctions in the heavy-quark limit, and represent their transverse 
momentum dependence explicitly. These results show that the dynamics within the 
two-particle Fock states is determined solely in terms of a single nonperturbative 
parameter A. 

Up to now, the transverse momentum distributions in the B mesons have been 
completely unknown, so that various models have been used in the literature. Some¬ 
times the dependence of the light-cone wavefunctions on the transverse separation 
is simply neglected, such as ^/>(w, —z^) = ip{uj, 0); clearly, this in general contradicts 
our results (14). Another frequently used models assume complete separation (fac¬ 
torization) between the longitudinal and transverse momentum-dependence in the 
wavefunctions, such as ip{ui,kT) = <j){u})T{kT) (see e.g. Refs.^^’^’^®). One typical 
example of such models is given by^ 

V’^^^(w,fcT) = A^w2(l - w)^exp X exp , (18) 

where N is the normalization constant and ujq = O.SGeV, K = 0.4GeV.'^ Our 
results (16) and (17) show that the dependence on transverse and longitudinal 
momenta is strongly correlated through the combination fc^/[u;(2A — w)], therefore 
the “factorization models” are not justified.'^ We further note that many models 
assume^^’^49 Gaussian distribution for the fer-dependence as in Eq. (18). These 
models show strong dumping at large \zt\ as ~ exp /2). In contrast to this 

asymptotic behavior, our wavefunctions (14) have slow-damping with oscillatory 
behavior as 4’±{u!, —z^) ~ cos(|2:t|\/w(2A — uj) — tt/A)/^J\ zT\■ 

To summarize, in this work we have derived a system of differential equations for 
the B meson light-cone wavefunctions including the transverse degrees of freedom, 
and obtained the exact solution within the valence Fock states. The differential 
equations are derived from the exact equations of motion of QCD in the heavy- 
quark limit. The heavy-quark symmetry plays an essential role as in the case 
of the light-cone limit. Heavy-quark spin symmetry reduces the number of 
independent wavefunctions drastically, so that the valence Fock-state components 
in the B mesons are described by only two light-cone wavefunctions. As a result, 
a system of four differential equations from the equations of motion becomes a 
“complete set” to determine these two wavefunctions, and this enables us to obtain 
the exact solution with full account of the fex-dependence. Also due to the power 
of heavy-quark symmetry, our final results are given in simple analytic formulae 
involving one single nonperturbative parameter A. Heavy-quark symmetry also 
guarantees that the solution in the present paper provides complete description of 
the light-cone valence Fock wavefunctions for the B* mesons and also for the D, 
D* mesons in the heavy-quark limit. 

■^In this model, the two independent wavefunctions of the B mesons are set equal to each 
other, . 

*The “non-factorization” in the light-cone wavefunctions for the light-mesons has been discussed 
in Ref. ^9, where the coupling between transverse and longitudinal momenta through the variable 
fc|,/[u(l — u)], with u the momentum fraction of the light quark, has been demonstrated. 
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Finally a comment is in order concerning the error induced by the Wandzura- 
Wilczek approximation. From the study of the B meson distribution amplitudes in 
the light-cone limit, there has been indication that, in the heavy-light quark systems, 
the higher Fock states could play important roles even in the leading twist level. 
This would suggest that the shape of the wavefunctions as function of momenta 
and their quantitative role in the phenomenological applications would be modified 
when including the higher Fock states. For example, inspecting the t ^ 0 limit of 
Eqs. (8) and (11) of Ref.^°, one immediately obtains the exact result for the first 
moment of as 


dujcfkT ^ t ) = 4 




dz^ 


- 3 + ^h) ) 


(19) 

where = 'ip± + denote the total wavefunctions including the higher Fock 

contributions and Xe and Xh denote the reduced matrix elements of rel¬ 

evant quark-antiquark-gluon operators in the notation of Ref.^°, representing the 
chromoelectric and chromomagnetic fields in the B meson rest frame, respectively. 
The first term |A^ in the RHS of Eq. (19) coincides with the moment of our solu¬ 
tion ipE, Eqs. (16) and (17), in the Wandzura-Wilczek approximation, while other 
terms | (X^ + come from the higher Fock contributions ■ The result 

(19), combined with an estimate X'^/A'^ = 0.36 ± 0.20, A|^/A^ = 0.60 ± 0.23 by 
QCD sum rules,suggests that the higher Fock contributions might considerably 
broaden the transverse momentum distribution. This point can be studied in a sys¬ 
tematic way, and more sophisticated wavefunctions will be discussed in detail 
in a separate publication.^^ However, qualitative features revealed in this paper, 
like non-factorization of longitudinal and transverse directions, “slow-damping” for 
transverse directions, etc., will be unaltered by the effects of multi-particle states, 
and helpful in elucidating QCD factorization theorems. 
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